Implications of a Running Spectral Index for Slow Roll Inflation 
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We analyze the weak {2a) evidence for a running spectral index seen in the three-year WMAP 
dataset and its implications for single field, slow roll inflation. We assume that the running is 
comparable to the central value found from the WMAP data analysis, and use the Hubble Slow Roll 
formalism to follow the evolution of the slow roll parameters. For all parameter choices consistent 
with a large, negative running, single field, slow roll infiation lasts less than 30 efolds after CMB 
scales leave the horizon. Thus, a definitive observation of a large negative running would imply that 
any inflationary phase requires multiple fields or the breakdown of slow roll. Alternatively, if single 
field, slow roll infiation is sources the primordial fiuctuations, we can expect the observed running 
to move much closer to zero as the CMB is measured more accurately at small angular scales. 



INTRODUCTION 

One of the most perplexing features of the 3-year 
dataset^ from the Wilkinson Microwave Anisotropy 
Probe [WMAP] 0, 0, is that it hints at a significant 
"running" in the scalar spectral index. This feature was 
seen, with roughly the same level of significance, in the 
WMAPI analysis [J, |a| ■ The evidence for a running in- 
dex persists when large scale structure data is added to 
the analysis ^ , but is diluted by the addition of Lyman- 
a forest data, both with WMAPI @ and WMAPII f^. 
Joint analyses require a melange of datasets to increase 
their coverage in fc-space and thus harbor the possibil- 
ity of systematic normalization uncertainties or a tension 
in relative normalization between datasets, which could 
manifest itself as a spurious running. This source of er- 
ror is eliminated when WMAPII is considered on its own, 
but only at the cost of reduced fc-coverage. 

We stress that WMAPII does not demand a non- 
vanishing running, as zero is ~ 2ct from the central value. 
Rather, we spell out the consequences of a large running 
- particularly with a negative sign - and the likely con- 
sequences for inflation. We find that a running similar to 
the WMAPII centroid rules out all simple models of infla- 
tion: those driven by a single, minimally coupled, scalar 
field whose evolution is well described by the slow roll 
formalism. Inflation, if it happens at all, would thus be 
non-minimal. Alternatively, future analyses of datasets 
with a larger fc-coverage will yield a running much closer 
to zero than the value extracted from WMAPII. 

We proceed by using the Hubble slow roll [HSR] hier- 
archy [a IE llO, ll3i ll^ J which connects astrophysical de- 
terminations of the primordial power spectra (both ten- 



sor and scalar) to the inflationary potential. In '12!] we 
showed how a Monte Carlo Markov Chain [MCMC] anal- 
ysis of cosmological data can constrain the HSR param- 
eters. These parameters are associated with flow equa- 
tions that determine their scale dependence. We use this 
system of equations to compute the number of efolds of 
inflation which occur after a mode with wavenumber fco 
leaves the horizon. Choosing ^ large enough to explain 
the central value found for the running in Spergel et al. 
guarantees that any single field model of slow roll infla- 
tion will run for less than 30 e-folds after CMB scales 
leave the horizon. As we explain below, one can increase 
the number of e-folds by adding a fourth slow roll pa- 
rameter. This term would be third order in the slow roll 
expansion,^ and the flrst three terms in the expansion 
would have roughly the same level of significance. Con- 
sequently, the slow roll expansion cannot be safely used 
to describe inflationary models with a substantial nega- 
tive running. 



THE SLOW ROLL APPROXIMATION 

As summarized in |l2| , the dynamics of single field in- 
flation can be written in the Hamilton- Jacobi form, where 
overdots correspond to time derivatives and primes de- 
note derivatives with respect to (j). The HSR parameters 
^Xh obey the infinite hierarchy of differential equations 
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1 We denote this WMAPII, and the one year dataset WMAPI. 



^ The usual counting is that e and tj are the lowest order terms, ^ 
second order, and the next term is thus third order. 
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The usual slow roll parameters are 77 = ^Xh and ^ = ^Xh- 
If we set the higher order terms to zero at some fiducial 
point, these differential equations ensure they vanish at 
all other times. The potential is given by 
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Liddle showed that the hierarchy can be solved exactly 
when truncated at order M IllL so 
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The Bi are specified by the initial values of the HSR 
parameters, 
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where the subscript refers to their values at the moment 
the fiducial mode fco leaves the horizon, and (j) = (j)o = 0. 
The number of e- folds, N is given by 

dN 47r H 
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Finally, (j) and k are related by 
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We now turn to the inflationary observables. 
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where C = 4(ln2 + 7) — 5, and we have introduced the 
customary notation a — dng/dlnk and r is tensor:scalar 
ratio. We retain all terms in a up to quadratic order in 
the slow roll parameters, anticipating that ^ may be as 
large as e or 77. 

COSMOLOGICAL CONSTRAINTS 

Consider the following very weak constraints 

0.9<n, <1.1 (11) 

r<l. (12) 
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FIG. 1: We show the regions of the (e,7;) plane excluded 
by the assumed bounds on the spectral parameters. As ^ 
increases, the portion of the plane for which a > shrinks. 



These bounds are broad enough to be independent of any 
assumptions about the presence or absence of a running 
index, and are outside the 99% confidence intervals de- 
rived from WMAPII and other datasets J3 • The running 
is given by equation 1101 - ^ is the only term linear in 
the slow roll parameters and it dominates the expression 
when I a I is large. In particular, for large negative a, we 
need ^ > 0. While e and 77 can become large for fixed ris^ 
they make a positive contribution to a thanks to the co- 
efficients on 677 Etnd e^ in equationEB FigureHshows the 
cuts these constraints put on the (e, r\) plane for ^ = 
and ^ = 0.01. Finally, we show the region where a > 0. 
Again we remind the reader that WMAPII does not 
require that a < and the purpose of this Letter is to 
explore the consequences of measuring a value of a near 
the WMAP centroid. This number depends on whether 
one allows for a contribution from primordial tensors. 
For the pure scalar case, the distribution peaks at a ~ 
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FIG. 2: We show the regions of the (e, 77) plane excluded by the assumed bounds on the spectral parameters for the same cuts 
displayed in Figure Q We include contours (for A'' = 15, 30, 45, 60) showing the number of e-foldings occuring after the fiducial 
mode fco ~ 0.002Mpc~^ leaves the horizon. For ^ = 0.03 the only regions of the (e, 77) plane with an acceptable value of n^ have 
N < 15. The running is a weak function of e and t; - the values of a above are taken from ^ alone. 



—0.05. Including tensors actually makes the central value 
larger, but this is not necessary for us to establish the 
conclusions we reach belcw. 

The minimal required amount of inflation is not ■well 
defined. If inflation is to happen before the electro^weak 
phase transition, the number of e- folds A^ must be greater 
than 30, and A^ ~ 55 for GUT scale inflation. These 
numbers are mildly dependent on the reheating mech- 
anism, but are sufficient for our purposes. Any model 
■with A'^ < 30 is unlikely to provide a ■workable explana- 
tion for the large scale homogeneity and isotropy of the 
universe. Because we can compute the running of the 
HSR parameters as a function of (j) (or, equivalently, fc), 
we can obtain the remaining number of e-folds for any 
choice of {eo, ryo, ^o}- The end of inflation is signified by 
the instant ■when e = 1, and we find N from equation ® 
alongside the fio-w equations, QJ and ^. The result of 



this calculation is displayed in Figure |21 for four values of 
^. Truncating the slo^w roll hierarchy at third order and 
choosing ^ consistent ■with a large negative running, we 
see that the fraction of parameter space that has A^ > 30 
and a < —0.02 is of measure zero. 



DISCUSSION 

In the above analysis we have sho'wn that a large, nega- 
tive of a, the running of the scalar spectral index, cannot 
be produced by a self-consistent single field infiationary 
model ■which is ■well described by the slo^w roll expansion. 
Specifically, the large running leads to an unacceptably 
small value of A^, the number of e-folds. The lo^w value 
of A^ can be ameliorated by adding a further slo^w roll 
parameter (^A) to the analysis ■which, if carefully cho- 



sen, can move ^(0) closer to zero as inflation continues. 
However, in this case the inflationary dynamics would be 
described by four numbers, {cq, tjq, ^q,^ Aq}, all of roughly 
equal importance. Such a model may be "slowly rolling" 
in the sense that e ^ 1, but the slow roll expansion is 
not trustworthy when its first three orders are of roughly 
equal weight. Nor is there is any guarantee it could be 
safely truncated after '^Aq. 

If the central value of a derived from the WMAPII 
analysis (with or without large scale structure informa- 
tion) is near the actual value, all models of inflation 
driven by a single, minimally coupled, slowly rolling 
scalar field are ruled out, and a large class of inflationary 
models would thus be falsified. There is no guarantee 
that inflation is well described by slow roll. Models with 
"features" in the potential yielding a local violation of the 
slow roll conditions have been considered in the past (e.g. 
[13, |lj, Il3)- Further, models of inflation with two or 
more interacting flelds have much more freedom to yield 
complicated spectra [ig, or one could turn to scenarios 
with two or more bursts of inflation [H [S E HE IHl . 
However, invoking any of these options would dramat- 
ically complicate the theoretical understanding of infia- 
tion. 

This analysis depends on the HSR equations, and their 
ability to incorporate constraints arising from the dura- 
tion of infiation alongside those from the perturbation 
spectra [l3|- The role of the ^ (or "jerk") term in the 
dynamics of the flow equations is discussed in detail in 
[23 . This analysis focuses on "attractors" in the slow roll 
parameter space, in the presence of non-zero high order 
terms in the slow roll hierarchy, rather than their end- 
state, which is our principal concern here. Likewise, in- 
flationary models generated via the flow equations [sl, |ia| 
or Monte Carlo Reconstruction [lOj can provide solutions 
with a large, negative running and 7V>55, but these rely 
on contributions from higher order terms. The accuracy 
of a specific class of slowly rolling models with signifi- 
cant running is probed using Monte Carlo Reconstruction 
techniques in [23j. Conversely, '24| argues that the sup- 
port for running in the WMAPI dataset comes from the 
multipoles with 30 < £ < 40. Finally, [l^ also point to 
the tension between a large running and providing the re- 
quired number of e-folds. However, our conclusion: that 
no single field model can simultaneously a) satisfy slow 
roll, b) produce a large negative running and c) provide 
a sufficient number of efolds, is new. 

Of course, by far the most prosaic resolution of this 
conundrum would be that the weak evidence for a signif- 
icant running in the WMAPII dataset is a statisical arti- 
fact that will evaporate as more data becomes available. 
Indeed, adding Lyman-a forest data is known to reduce 
the need for a negative value of a [a,l2l- In particular, we 
anticipate that the Planck mission will significantly en- 
hance our understanding of a by providing high quality 
measurements of the fundamental power spectrum over a 



larger wavelength range than WMAP, without the need 
to combine several heterogeneous datasets that carry the 
risk of systematic biases in their relative normalization. 
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